Exact sampling and entanglement-free resources for measurement-based quantum computation 
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We identify a family of quantum computations whose output is likely hard to exactly simulate (sample) clas- 
sically. Specifically, if the output distributions can be efficiently and exactly simulated classically then the 
polynomial hierarchy from computational complexity theory collapses at its third level. We demonstrate that 
these circuit familes can be implemented in MBQC using resources states with zero entanglement and zero 
quantum discord, that is, states which are essentially classical probability distributions. This work shows that 
families of states exist which violate no Bell inequality, but which are nevertheless unlikely to be efficiently re- 
alizable in any theory respecting classical notions of computational complexity — such correlations thus possess 
a complexity theoretic nonclassicality. 



There is a strong belief that quantum computers can ef- 
ficiently perform certain tasks that cannot be performed ef- 
ficiently on a classical computer, such as integer factoriza- 
tion [ 1 1 . One of the central problems of quantum informa- 
tion theory is to determine the origin of the apparent quantum 
speed-up. 

There has been much interest in this question since the early 
years of quantum computation. For example, Jozsa and Lin- 
den showed that increasing multi-partite entanglement is nec- 
essary for a pure state quantum computer to achieve an expo- 
nential speed-up 0. In contrast, it has long been suggested 
that mixed state computation may not require unbounded en- 
tanglement for a quantum speed-up. It is speculated that quan- 
tum discord may play an important role instead, for example, 
in the "one-clean qubit" model [3, 4J. The subtlety of this 
question is illustrated by the Gottesman-Knill theorem, which 
says that Clifford circuits, which can produce highly entan- 
gled states, can nevertheless be efficiently simulated on a clas- 
sical computer [5 1. Furthermore, Van den Nest has shown that 
a quantum computer can retain its full computational power 
even when certain measures of entanglement are restricted to 
be vanishingly small [6|. 

Investigating this question is made all the more difficult by 
the lack of a rigorous proof that quantum computing does in- 
deed have a computational advantage over classical compu- 
tation. For example, there is no proof that an efficient clas- 
sical factoring algorithm does not exist. Nevertheless, rigor- 
ous computational complexity theoretic evidence can be given 
that if a classical computer could exactly simulate a quan- 
tum computer |7| — and by "simulate" we mean reproduce the 
same input-output statistics and thus sample from the same 
probability distribution — there would be very significant con- 
sequences for computational complexity theory which violate 
most experts' intuitions. There have been a number of notable 
examples of such evidence, in particular by Bremner, Josza 
and Shepherd (BJS) [8| and Aaronson and Arkipov (AA) J5). 
Surprisingly, the computations they considered are not univer- 
sal models, but much simpler models such as Instantaneous 
Quantum Polytime (IQP) circuits, and linear optical boson 
sampling. The evidence provided by BJS and AA rests upon a 



"generic, foundational" assumption of computer science [9|, 
namely that the levels of the polynomial hierarchy (PH) ifTOll 
are distinct, in particular that the third level of the hierarchy is 
strictly smaller than the highest level. 

Hypothesis 1. The third level of the Polynomial Hierarchy is 
strictly smaller than at least one other level in the hierarchy. 



This hypothesis is similar in spirit to the common assump- 
tion P 7^ NP. BJS and AA showed that if IQP or boson sam- 
pling circuits can be efficiently and exactly simulated with 
classical computers then the PH collapses to its third level 
(meaning that all levels higher than the third level are identi- 
cal to the third level), violating Hypothesis 1 . 

The simplicity of AA and BJS's models raise the hope that 
a quantum experiment might soon demonstrate a computation 
beyond the power of conventional computers. Early experi- 
mental demonstrations of boson sampling ifTTl have been in 
agreement with quantum mechanical predictions. 

In this paper, we provide further evidence of the subtle na- 
ture of quantum speed-up by showing that there exist com- 
putations which a quantum computer can perform more effi- 
ciently than a classical computer without using either entan- 
glement or discord. 

We do this by modifying the argument presented by BJS 
to derive a smaller set of uniform families of circuits (IQP*) 
and show that a similar hardness of simulation theorem holds 
for this family. We then use this to derive two main results, 
which both hold if Hypothesis 1 is true. Firstly, we show that 
families of separable and discord-free states exist for which 
efficient and exact sampling from measurements in the com- 
putational basis is impossible. Secondly, we show that there 
exist computations which cannot be efficiently and exactly 
performed on a classical computer, but can be performed in 
standard measurement-based quantum computation (MBQC) 
using resource states with zero entanglement and zero discord, 
and consequently with correlations that do not violate any Bell 
inequalities. 

These results show that there exist classical probability dis- 
tributions which can be efficiently generated by simple quan- 
tum computers which nevertheless are unlikely to be effi- 
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FIG. 1: a) Standard form of an IQP circuit, where each gate is diag- 
onal in the Pauli-X basis and (3:1,0:2) is the two bit input string. All 
measurements are in the computational (Pauli-Z) basis. The boxed 
gates give the unitary D. b) a MBQC implementation of the circuit 
in a), where each circle represents a qubit prepared in the state |+) 
and edges between circles represent the application of a controlled- .Z 
gate. The contents of each circle represent the basis in which the cor- 
responding qubit is measured, where X represents the Pauli-X basis, 
and 6j represents the basis Ux{—0j)ZUx{0f), where Ux(8j) is a 
rotation by 9j about the Pauli-X axis. The angles Oj are in one-to- 
one correspondence with the 9 Z in the representation of Dj in Eq. 
(TV All of these measurements can be implemented simultaneously 
(non-adaptively) in MBQC. 



ciently realizable by a classical device. Furthermore, correla- 
tions encoding such distributions are computational resources 
(e.g. for MBQC), allowing one with access to them to effi- 
ciently perform apparently nonclassical computations. This 
provides an example of a form of nonclassicality for probabil- 
ity distributions derived from complexity theoretic arguments 
rather than via hidden variable models. 

We begin by defining IQP (Instantaneous Quantum Poly- 
time) circuits fl2l . introduced by BJS, which play a central 
role in this paper. 

Definition 1. An IQP circuit with classical input bit string x 
of size n acting on q > n qubits consists of: 

1. a quantum register prepared in the input state 

\x)\0f q ~ n ;<md 

2. the application of a unitary operator U to the register, 
where U is diagonal with respect to the eigenbasis of 
Pauli-X operators. 

We denote the output of this computation, obtained via com- 
putational basis measurements on every qubit, by the g-bit 
string m, whose jth element rrij e {0, 1} is the outcome of 
the computational basis measurement on the jth qubit. 

An example of an IQP circuit is illustrated in Fig.fTk. Such 
circuits are called instantaneous because D can be decom- 
posed into a product of commuting gates, which can thus be 
applied in any order (or simultaneously) [12|. 

In studying the computational power of a family of circuits, 
one must ensure that unreasonable computational power is not 
hidden in the description of the circuits themselves. This can 
be ensured via a uniformity condition which ensures that a de- 
scription of each circuit in the family can be efficiently gener- 
ated. While BJS introduce a uniformity condition in Ref. [8|, 
it is both convenient and necessary for our analysis to adopt 
a different one here, and we denote the set of uniform circuit 
families obtained under this condition by IQP* . 



Definition 2. An IQP* circuit family is a family of IQP cir- 
cuits, with input x and input size n = \x\, followed by com- 
putational basis measurements on every qubit, such that the 
number of qubits q is polynomial in n, and where the unitary 
operator U n (which has an explicit n -dependence) is a poly (n) 
product of gates of the form 



D{9 Zl z) 



,x\. 



(1) 



where angle 9 Z <E (0, 2tt] has a description polynomial-size in 
n, z is a g-bit string, and we introduce the notation X[z] — 
■ X z f where Zj is the jth bit of z. E.g. D(0 W1 , 101) = 
exp[i0ioi(X(g) 1 ®X)]. 

In other words, the description of every member of an IQP* 
circuit family is a polynomial list L n of g-bit strings z and 
corresponding angles 8 Z . This list then defines the circuit for 
input size n. If we adopt the notational shorthand that 6 Z = 
for all bit strings z not in L n , then the unitary transformation 
for the circuit U n is given by 
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Operators of this form have a useful symmetry in their matrix 
elements, which we will exploit below, namely, 



(w®y\D n \w) = (y\D n \0) 



(3) 



for all bit strings w and y, with |0) = |0) q and where © rep- 
resents a bit-wise sum modulo 2. 

This is different from the uniformity condition introduced 
by BJS [8] in some significant ways. Most importantly, the 
circuit families in [8| allow the circuit to depend on the input 
string x, rather than the length of x. This makes BJS's set of 
uniform families much larger than IQP*. Indeed BJS's set of 
uniform circuit families include all polynomial-sized classical 
circuit families, but IQP* circuit families cannot even achieve 
a deterministic Toffoli gate. Consequently, BJS's hardness of 
simulation theorem does not immediately apply to IQP* cir- 
cuit families. We prove here that even this much smaller set 
of uniform circuit families are unlikely to be efficiently and 
exactly simulatable on a classical device. 

Lemma 1. The output probability distributions generated by 
IQP* circuit families cannot be efficiently and exactly simu- 
lated on a classical computing device unless Hypothesis 1 is 
false. 

The full proof of Lemmafllis presented in Appendix[B| The 
technical definition of a classical computing device is pro- 
vided in Appendix [A] along with other useful notions from 
computational complexity theory. To summarize the proof for 
readers familiar with Ref. |8), under postselection of mea- 
surement outcomes of a subset of qubits, the families of IQP* 
circuits can be mapped to general quantum circuits satisfy- 
ing the standard uniformity condition for the complexity class 
BQP. This allows us to utilize a similar proof technique to BJS 
in Ref. (8). 

As in Ref. (8), Lemma [T] may be generalized to include 
multiplicative error up to a certain factor on the individual 



probabilities. However, since such errors are not physically 
relevant, we omit the generalization. Lemma[T]presents strong 
evidence that IQP* circuit families may not be efficiently and 
exactly simulated on a classical computer. 

We now use Lemma [T] to show that computational basis 
measurements on the following set of states can cannot be ef- 
ficiently and exactly simulated on a classical computer. 

Definition 3. An IQP* zero-input state family is the set of 
quantum states created by an IQP* circuit family, where the 
input is set to the all zeros string ... 0. 

Theorem 1. The statistics of computational basis measure- 
ments on IQP* zero-input state families cannot be efficiently 
and exactly simulated on a classical computer unless Hypoth- 
esis 1 is false. 

Proof We use a special property of IQP circuits, namely that 
the output statistics of an IQP circuit with input x, defined 
according to Definition 1, may be realized by the same IQP 
circuit acting on the rc-bit all-zeros string by performing some 
simple extra post-processing of the output bits of the measure- 
ments. 

Observe that the probability that the measurement output 
string is a bit string m given input x is 

Prob.(77i|x) = \(m\D n \x}\ 2 , 

= \(m®x\D n \0)\ 2 , 

= Prob.(m©x|0), (4) 

where a; is a; appended by q — n zeros and we obtained the 
second line from Eq. (BJ. 

Thus identical output statistics to an IQP circuit given input 
x can be obtained via the same circuit with input and post- 
processing of the output bit string m to string m © x. This 
post-processing comprises at most n bit-flips and can be (triv- 
ially) efficiently performed on a classical device. From this, 
the theorem follows directly from Lemma 1 . □ 

Consideration of dephasing with respect to the computa- 
tional basis provides an immediate corollary. 

Corollary 1. There exist uniform families of separable 
discord-free quantum states (IQP* state families completely 
dephased in the computational basis) that when measured in 
the computational basis produce statistics that cannot be effi- 
ciently and exactly simulated on a classical computer unless 
Hypothesis 1 is false. 

Depolarized states can be considered equivalent to classi- 
cal probability distributions, so Corollary [TJimplies that there 
exist classical probability distributions that can be efficiently 
generated via quantum means which are likely to be impossi- 
ble to efficiently generate classically. 

We now show that entanglement-free and discord-free 
states can be a computational resource by considering the im- 
plications of Lemma[T]for measurement-based quantum com- 
putation (MBQC) O [14). In all MBQC schemes known 
to date, the resource states are entangled. Furthermore, in 
models of MBQC where the classical side-processing only in- 
volves computing bit string parities alone, (e.g. cluster state 



and graph state approaches) it has been shown that the viola- 
tion of Bell inequalities (and hence entanglement) is necessary 
for universal quantum (indeed universal classical) computa- 
tion lfT5ffl"8l . Armed with Lemma 1, however, we can show 
that MBQC can utilize resource states with neither entangle- 
ment nor quantum discord to achieve a computation which no 
classical device can exactly and efficiently simulate. 

Theorem 2. IQP* circuit families may be implemented in 
MBQC using resource states with zero entanglement and zero 
discord. 

The proof of this theorem, presented in full in Appendix 
[C] rests on the following two facts. Firstly, any IQP circuit 
can be implemented in MBQC using solely non-adaptive mea- 
surements |[T2l[T9ll20l ; an example of an implementation is 
given in Fig. [Th Secondly, given a resource state \i/j) and a 
fixed set of local projective measurements on all the qubits of 
\ip), we can apply a fully dephasing channel in the respective 
measurement eigenbasis for each qubit without changing the 
measurements statistics. Any state p which is the output of a 
completely dephasing channel on every qubit is separable and 
discord-free [3|. 

Theorem|2]immediately leads to the following corollary. 

Corollary 2. IQP* circuit families may be implemented in 
MBQC using correlations which do not violate any Bell in- 
equalities. 

Proof Fully dephased states, diagonal with respect to a local 
basis, provide an immediate local hidden variable model that 
can simulate the MBQC implementation of an IQP* circuit. 

□ 

Discussion. Throughout this paper we have used the as- 
sumption of Hypothesis 1, that the Polynomial Hierarchy of 
complexity classes does not collapse to its third level. Given 
this assumption we have shown that there exist families of sep- 
arable discord-free quantum states on which computational 
basis measurements cannot be efficiently and exactly simu- 
lated on a classical computer. Furthermore, closely-related 
families of such states can be resources for performing IQP* 
computations in MBQC. 

The statistics of computational basis measurement out- 
comes define a family of classical probability distributions. 
Thus these results imply that certain classical probability 
distributions that are unlikely to be exactly and efficiently 
achievable with classical devices (although they have an effi- 
cient quantum implementation) may themselves be resources 
for nonclassical computation, analogous to the entangled re- 
sources states for universal MBQC. These probability distri- 
butions violate no Bell inequality, but have the nonclassical 
property that they are unlikely to be realizable in any the- 
ory respecting classical notions of computational complexity. 
This can be considered a form of nonclassicality for correla- 
tions based on complexity theoretic restrictions alone. 

The output probabilities for random linear optical networks 
(boson sampling) |9| share this property, and, based upon the 
simplicity of the quantum processes that generate such prob- 
abilities, we conjecture that this complexity-theoretic form of 
nonclassicality is common in interacting quantum systems. 



Theorem 2 illustrates that discord-free and entanglement- 
free states, i.e. (in essence) classical probability distribu- 
tions, are resources for input-dependent families of nonclas- 
sical computations. This shows that these distributions have a 
resource character not known to be exhibited by boson sam- 
pling distributions. In the context of MBQC, they provide 
computational resources for nonclassical computation. Thus, 
not only are the distributions hard to simulate, they are re- 
sources for a nonclassical computational task. To achieve this 
in boson sampling, a uniform family of computations with an 
input would need to be formulated. Whether this is possible 
we leave as an open question — the links identified between 
boson sampling and IQP |9| are promising in this regard. 

It is important to consider why these distributions cannot be 
efficiently classically generated. The complexity of sampling 
from general probability distributions via a classical computer 
given a source of uniform random bits was first studied in de- 
tail by Knuth and Yao in 1976 [21]. They proposed a sampling 
algorithm, for any probability distribution over n-bit strings, 
with an average run-time less than n + 2. However, to set up 
their algorithm in general requires exponential computational 
resources and has an unbounded worst-case runtime, meaning 
that for certain probability distributions, a finite exact simula- 
tion via their algorithm is impossible. 

It is also natural to ask how general these probability distri- 
butions are — can all probability distributions be achieved via 
an (appropriately large) IQP circuit? We show in appendix [D] 
that this is not the case. For the special case of a two-qubit 
IQP circuit, the range of distributions achieved is broad, but 
does not cover all two-bit distributions. 

We should emphasize that these results do not imply that 
nonclassical computations can be performed using quantum 
circuits that do not create discord or entanglement at any time- 
step. Indeed, the results of Eastin [22 1 show that unitary 
discord-free computations (consisting of one- and two-qubit 
gates on a restricted family of input states) are classically ef- 
ficiently simulatable. This implies that the resource states in 
Theorem [2] cannot be efficiently generated via a unitary cir- 
cuit and then applying a dephasing channel without creating 
entanglement in an intermediate step. Methods for the genera- 



tion of such states with a minimum amount of nonseparability 
remains a topic for future work. 

One important limitation of our method, similar to that of 
Ref. (HI, is that our theorems only consider exact simula- 
tion. Our theorems would have an even broader impact if 
they could be shown to be robust to error, as this would allow 
for approximate simulation. While the limited approximation 
model (multiplicative errors on individual probabilities) con- 
sidered in Ref. [8| can be immediately applied here, it does 
not provide an analysis of the most physically reasonable error 
model, namely, additive error on the whole probability distri- 
bution. We note that AA (9) have made further progress in this 
direction in the context of boson sampling, and in further work 
we shall aim to provide evidence of hardness for reasonable 
approximate simulations of the quantum processes considered 
here. The mapping of this problem to a solely classical simu- 
lation problem may aid in this research. 

This work shows that even very simple physical quantum 
models can retain a classical hardness of simulation, motivat- 
ing further investigation into such models. If such a model 
arose in a natural physical system which we wanted to simu- 
late, this would provide a powerful application of these sim- 
pler quantum computers. 

The relationship between entanglement and quantum com- 
putational speed-up has been debated since the early days of 
quantum computing theory. These results emphasize that, 
when considering computation with mixed states, we can- 
not make assumptions on the necessity for either entangle- 
ment nor quantum discord. The relationship between quan- 
tum computational power and correlations appears even more 
subtle than previously thought. 
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Appendix A: Some computational complexity definitions 

It is important that we define what we mean by "classical 
computing devices". By this expression we include general 
probabilistic processes since quantum circuits will, in general, 
not give deterministic outcomes and we want to simulate the 
statistics of quantum circuits. We follow [8 1 in making the 
following definition (again where Boolean circuits are defined 
in (TO)): 

Definition 4. A classical computing device takes an input 
x £ {0, 1}™ of size n and produces a bit string y e {0, 1} S of 
size s by performing the following operations: 

1. Flip n fair coins to produce a random bit string z £ 
{0, l} ri of size r\ where we make the assignment 
"heads" to and "tails" to 1; 

2. Prepare the bit string x' = (x, 0, z) where is the bit 
string of size r 2 with all elements equal to 0; and 

3. Apply a Boolean circuit B n that takes bit string x' of 
size n + r\ + r 2 as an input and outputs the bit string 
y G {0, 1} S of size s, where the description of B n is 
generated in poly(n) by a classical Turing Machine. 

All variables n, s, r\, r 2 are positive integers. If the com- 
putational resources of this device are polynomial in n then 

s = poly(n), r\ = poly(n), and r 2 = poly(n). 

There is a classical computational complexity class asso- 
ciated with these devices if the computational resources they 
use are polynomial in n, and this is the class of decision prob- 
lems known as BPP. Since these classical computing devices 
involve nondeterministic processes, the correct answer to a de- 
cision problem may not be obtained deterministically; there is 
a probability of giving the wrong answer. A problem is in BPP 
if this error probability is less than or equal to some constant 
c < 1/2. A related and potentially larger complexity class 
of probabilistic classical computations is PP where the error 
probability is bounded from above by 1/2 but may not be a 
constant, i.e. can depend on the input size n. This last class 
will also be important in later discussion. 

Similar to BJS [8| and AA |9|, our proof makes use of com- 
plexity classes defined using post-selection. Post-selection 
cannot be achieved deterministically in a physical realization, 
but is an extremely useful technical concept. Post-selection 
is the act of demanding that the outcome of a quantum mea- 
surement is a fixed value, and that the state of the system then 
evolves via a fixed projector and a renormalization. Aaron- 
son introduced the class PostBQP, which, loosely speak- 
ing, represents the problems solvable on a quantum computer 
given polynomial resources if we were given the extra power 
of post-selecting measurement outcomes. BJS introduced the 
complexity class PostIQP, which applies a similar treatment 
to their uniform IQP circuits. 

The formal definitions of the classes PostBQP, PostIQP 
and PostIQP* are as follows. 

Definition 5. A BQP circuit family {C„ : n e N} is a set of 
quantum circuits such that for each input bit string x of size 



n £ N, C n is a quantum circuit acting on q = poly(n) qubits 
(initiated in the state |x)|0)® 9_n ) and is a sequence of gates 
chosen from the universal gate set {CZ, H, Z, P} and has a 
description generated in poly(n) time by a classical Turing 
Machine. 

Here CZ is the controlled-./? gate, H is the Hadamard gate, 
Z is the Pauli-Z gate and P = e l i z . We now define the afore- 
mentioned complexity classes under post-selection utilizing 
the definitions above and following the work of Aaronson in 



Definition 6. A language L C {0, 1}* is in PostBQP iff 
there exists a BQP circuit family {C n \n £ N} such that for all 
inputs x £ {0, 1}": 

1. after C n is applied to the state |x)|0)|0)...|0), the last 
qubit has a nonzero probability of being measured in 
the computational basis as |1); 

2. if x £ L and the last qubit is in the state |1), the second- 
to-last qubit when measured in the computational basis 
is |1) with probability > 2/3; 

3. if x $ L and the last qubit is in the state 1 1) , the second- 
to-last qubit when measured in the computational basis 
is |1) with probability < 1/3; 

where the last qubit is the bottom qubit line in a quantum 
circuit diagram such as in Fig. [T] 

To define PostIQP, BJS use a near identical definition, 
replacing "a BQP circuit family" with their definition of "a 
uniform IQP circuit family" in Definition [6] Similarly, we 
define PostIQP*, by replacing, "a BQP circuit family" in 
this definition by "an IQP* circuit family". 



Appendix B: Proof of Lemma[T] 

In this section we provide a proof of Lemma[T] This section 
utilizes some of the technical concepts defined in Appendix [A] 
To prove LemmafTlwe use the arguments presented in (8). 

One of the key lemmas which underpin BJS's main result 
(see Corollary 3.3 in [8]) is the complexity class equation: 

Lemma 2. PostIQP = PostBQP = PP 

where PostBQP, PostBQP and PP are defined in Appendix 
[Aland in [8 1. The right-hand equation is due to Aaronson ll23l . 
To prove Lemma[T] we need to show that: 

Lemma 3. PostIQP*=PostBQP=PP. 

Lemma 1 then follows directly from all other steps of the 
proof of Corollary 3.3 in (8). We shall not reproduce those 
steps of the proof here. 

In order to prove PostIQP* = PostBQP, it is necessary 
to prove that PostIQP* C PostBQP and that PostBQP C 
PostIQP*. It is clear that the former is true. To prove the 
latter, recall that any BQP circuit can be expressed in the for- 
malism of MBQC [14|. The physical realization of any BQP 
circuit family in MBQC comprises the generation of a graph 



state of sufficient size and appropriate structure, followed by 
adaptive single-qubit measurements in the bases X, Y and 
{X ± Y)/y/2. 
Graph state generation comprises: 

1. preparation of a set of qubits in state |+) = (|0) + 
|l))/\/2;and 

2. implementation of CZ gates between certain qubits. 

The measurements can comprise rotations about the Z axis 
by any of the possible angles 0, ir/4, — 7r/4 or tt/2 (with re- 
spect to the Z basis) followed by a measurement in the X 
basis. 

The sequence of measurements corresponds to the chosen 
BQP circuit. Notice that from the definition of a BQP circuit 
there is a dependence on the value of x in the initial state and 
the gates in the circuit are generated as a function of n and not 
x. This dependence on x is introduced by having measure- 
ments of the following form (via a simple application of gate 
identities in [1141 ): there are n qubits that each correspond to 
each element Xj of x, and if Xj — we make the measurement 
in the X basis, and if Xj = 1 then we implement a n -rotation 
about the Z axis and measure in the X basis. Since the rota- 
tion commutes with the CZ gates, it can be incorporated into 
the preparation of that qubit by preparing |+) if Xj = and 
| — ) otherwise. For the rest of the qubits in the graph state, the 
rotation about the Z axis and measurement in the X basis is 
only defined by size n, the size of x. 

If we could post-select on measurement outcomes 
in MBQC then we remove the need for adaptive 
measurements — computations are accepted only if cer- 
tain outcomes occur. We now can show that computations 
in MBQC that no longer have adaptivity are instances of 
IQP* circuits. This can be easily shown if one writes out 
the computation in MBQC as a quantum circuit, and for all 
qubits after preparation of states |+) and |— ) and before mea- 
surements in the X basis two Hadamard gates are applied, 
yielding exactly the same circuit. The action of a Hadamard 
on preparation will map -75 (|0) + (— l) Xi |l)) — » \Xj). In a 
computation in MBQC all unitaries prior to measurement are 
diagonal in the Z basis, so the action of a Hadamard on every 
qubit prior to and after these unitaries maps these unitaries 
diagonal in the Z basis to unitaries that are diagonal in the X 
basis; those unitaries that appear in IQP* circuits. Finally the 
action of a Hadamard prior to measurement in the X basis 
results in a measurement in the computational basis. In Fig. 
[2] we give an example of this equivalence between MBQC 
circuits without adaptivity and IQP circuits. 

Therefore, under the action of post-selection MBQC com- 
putations are instances of IQP* circuits. To complete the proof 
of Lemma [T] we insert Lemma [3] into the proofs of Theorem 
3.2 in [8 1 and Corollary 3.3 in [8 1. 



Appendix C: Proof of Theorem 2 

The proof of this theorem can be broken down into three 
steps. First we show that any IQP circuit can be implemented 



1+) 




h\ — \rA\ 



FIG. 2: The top circuit where Dz is a unitary that is diagonal in the 
Pauli-Z basis is equal to the bottom circuit where Dx is a unitary 
that is diagonal in the Pauli-X basis. The top circuit is of the form 
of a circuit in MBQC but without adaptivity. Since H ■ H — 1 
we immediately see that the bottom circuit is an example of an IQP 
circuit. 



in MBQC without adaptive measurements |[T9l |20l , that is, 
where all measurement bases are fixed at the start of the com- 
putation. Next, we use some well-known properties of mea- 
surements and states: if any qubit on which observable A is 
to be measured undergoes a quantum channel where all Kraus 
operators commute with A, then both statistics and the post- 
measurement state are identical to the case where the channel 
was not applied. Finally, we show that a state undergoing a 
local dephasing channel on every qubit is both separable and 
discord-free. 

Lemma 4. Any IQP* circuit can be implemented in MBQC 
deterministically with non-adaptive, fixed-basis measure- 
ments. 

Proof. Recall from the proof of Theoremfllwe do not need to 
consider any dependence on the input x in the quantum circuit 
but only in post-processing of measurement data. Therefore, 
we only consider IQP* circuits where the input is the register 
|0)® 9 . As shown directly in Appendix [b] and demonstrated 
in Fig. [2] an IQP circuit can consist of the preparation of q 
qubits in the state | +) followed by a unitary that is diagonal 
in the Pauli-Z basis, and a measurement in the Pauli-X basis. 
In MBQC the preparation and unitary that is diagonal in the 
Pauli-Z basis can be implemented by: 1) preparing a graph 
state with r qubits for r > q and; 2) making measurements 
on r — q of those qubits. The remaining q qubits correspond 
to those qubits in the quantum circuit that now (as a result 
of the measurements) have undergone the necessary unitary 
evolution. 

To be more exact, Browne and Briegel 1 20 1 showed that any 
unitary gate of the form 



D z (6 z ,z) 



(CI) 



where Z[z] = (££) Z Zj and Zj is the jth bit of z, can be 
achieved in MBQC with simultaneous measurements, and 
with byproduct operators which are a tensor product of Z and 



they can be applied at the end of the computation and mea- 
surements do not need to be adaptive, that is, are not influ- 
enced by prior measurement outcomes. Therefore, the r — q 
measurements implement all sequences of unitaries CI in an 
IQP circuit and can be implemented simultaneously. 

All that remains to do is measure the remaining q qubits 
in the X basis to implement the IQP circuit. As Browne and 
Briegel demonstrate in |20|, these measurements can also be 
performed at the same time as the measurements on the r — q 
qubits that implement the desired unitary. After some simple 
classical post-processing of measurement data, an IQP circuit 
is implemented deterministically. Also, the choices of mea- 
surement and graph state prepared only depend on n, the size 
of the input x, and not x itself. □ 



Lemma 5. Given a state p, an observable O, and a state S(p) 
where S(-) is a superoperator whose Kraus operators all com- 
mute with O; then the probabilities of outcomes of the mea- 
surement of O on p and S(p) are equal. In addition, if p' is 
the post-measurement state given p, then S(p') is the post- 
measurement state given S(p). 

Proof. This is a well-known fact which arises from the Kraus 
operators and eigenprojectors of O being diagonal with re- 
spect to the same basis. □ 



Now consider the MBQC implementation of an IQP circuit 
introduced above. Let \ip) be the graph state and let Oj be the 
observables representing the projective measurement on the 
jth qubit. Every qubit in the state is measured. The channel 
S(p) with Kraus operators (1/^2)1 and (l/y/2)Oj, is called 
the dephasing channel where Oj is a qubit observable. Ac- 
cording to Lemma B] applying the dephasing channel to every 
jth qubit of \ip) will generate a mixed state p which, in terms 
of the the output statistics of measurements Oj, has identical 
properties to \t/j). Recall that this MBQC implementation al- 
ways has exactly the same measurements for every instance 
of x of size n. Hence, we can replace \ip) by p in the MBQC 
implementation of the IQP circuit. 

The final step in the proof of Theorem [2] is to demonstrate 
the (known) fact that p is separable and discord-free. As 
shown, e.g., in [3|, any discord-free n-qubit state p can be 
written 



Po 



= J2P°V U °\V)(V\ U 1 



(C2) 



1. An IQP circuit is then built of gates of the form CI Since 



the byproduct operators commute with the logical gates CI 



where <7 = ®" =1 Uq for Uq being a local unitary acting 
on the jth qubit, y is an n-bit string and the state \y) = 
\yi) <S> 1 2/2) <8> • • • \y n ) is a tensor product of computational ba- 
sis states. A few lines of algebra suffice to show that the state 
p generated by the action of the local dephasing channels on 
\ijj) is of the form of equation ( |C2| >. Therefore, for any IQP 
circuit, we can find a state p which is both entanglement- and 
discord-free, and a resource for the MBQC implementation of 
the circuit, completing the proof of Theorem [2] 



Appendix D: Two-qubit IQP Circuits and Classical Probability 
Distributions 

In this section we investigate, as a specific example, the 
conditional probability distributions generated by two-qubit 
IQP circuits. We show that while a broad range of probabil- 
ity distributions are possible with such circuits, they cannot 
generate all probability distributions over two classical bits. 
We are now in the framework of IQP circuits where the input 
of the circuit is the state |00) and we allow an arbitrary IQP 
unitary and measure in the computational basis. There is no 
dependence in the circuit on an input bit string x. 

Recall that the unitary in an IQP circuit can be described by 
a combination of bit strings to indicate which qubits are being 
acted on and the phase in the action on these qubits. These 
pairs are (z,6 z ) C ZJx (0, 2n] where n is now the num- 
ber of qubits in the circuit. The question is now whether this 
most general IQP circuit can produce a probability distribu- 
tion p(m) :— Prob.(m) over bit strings m G {0, l} n . We now 
demonstrate this with the most general IQP circuit for two 
qubits, thus corresponding to the generation of two classical 
bits (generally probabilistically). 

The most general IQP circuit for two qubits is composed of 
three gates Dg z for all bit strings z except (0,0). For clarity, 
we denote the three corresponding angles as 6\ = On, 82 = 
6*io and 6*3 = 6 J oi- The action of this circuit on input |00) is 




(0,0,0) 



u.o.or 



FIG. 3: Space of probability distributions over two bits including 
the simplex, the most general space of probability distributions. The 
light red points are statistics produced by random IQP circuits. The 
extreme points of the simplex are labeled by row vectors P as de- 
scribed in the text. The plane intersecting extreme points (0, 0, 1), 
(0, 0, 0) and the third point (1/2, 1/2, 0) demonstrates that the IQP 
circuits do not produce distributions that are on the face of the sim- 
plex. Also, as can be seen, IQP circuits cannot produce all statistics 
in the interior of the simplex. 



D|00) = Dg a D 02 D ei \OO) 

= [Q00I1I2 + «oiX 2 + a 1Q X 1 + 011X1X2] 1 00) 
= ooo|00) + a i|01) +aio|10) + ou|ll), 

where for notational convenience we have suppressed tensor 
products and Xj on its own is a Pauli-X acting on the jth 
qubit with the identity acting on the other. The amplitudes 
are: 



a-io 

an 



cos 
icos 
icos 
i sin 



#1) cos(6>2)cos(#3) — isin(#i) sin(#2)sin(#3), 

:os(9i) cos(# 2 ) sin(# 3 ) — sin(0i) sin(# 2 ) cos(6> 3 ), 
:os(9i) 8111(6*2) cos(^3) — sin(^i) cos^) sin(6>3), 
>in(6>i) cos(#2) 003(6*3) — 008(6*1) sin(6 2 ) sin(6 ) 3 ). 

When a computational basis measurement is made then 
\& mi m 2 1 2 is the probability of obtaining the classical bit string 
m = (mi,m 2 ). 

We now focus on probability distributions where p(0, 0) = 
0, thus implying that cos(9j) = sin(6*fc) = for j ^ k. If 
j = 1, then if k = 2 this necessarily means p(0, 1) = 
or if k = 3 then necessarily p(l, 0) = 0. Again by similar 
observations if j — 2, for k = 1 then p(0, 1) = and for 
k = 3 we necessarily have p(l, 1) = 0. Finally for j = 3, 
then necessarily p(l, 0) = and p(l, 1) = for k = 1 and 
k = 2 respectively. We could have started by setting any other 
p(m) to be zero and obtain the same observations. Therefore, 
IQP circuits cannot produce probability distributions where 
only one of the bit strings does not occur. That is, if one bit 
string never occurs then necessarily at least one other bit string 
cannot occur. 



If we form the simplex of probability distributions over two 
bits, we have a three-dimensional space of vectors P G M 3 
with elements p(m) for all m except m = (1, 1). All possi- 
ble probability distributions (including those produced by IQP 
circuits) can be written as a convex combination of the ex- 
treme points of this simplex, where each of the four extreme 
points corresponds to deterministic ally outputting a bit string 
m, i.e. p(m) = 1 for one m. While IQP circuits can attain 
these extreme points through appropriate choices of phases, 
they cannot produce probability distributions on the faces of 
this simplex; this would involve only one output m never oc- 
curring, which we have shown is not possible. In Fig. [3] 
we show statistics produced by randomly chosen two-qubit 
IQP circuits to point out how these circuits cannot produce 
all statistics for two bits. Since IQP circuits can achieve arbi- 
trary product distributions, these no-go distributions must all 
be nontrivially correlated. 

If we write the vectors P as row vectors 
(p(0, 0),p(0, l),p(l, 0)), then an example of a distribution 
that two-qubit IQP circuits cannot achieve is (1/4, 1/4, 1/2). 
A simple classical machine could produce this distribution by 
flipping two coins and translating "heads" as and "tails" as 
1: flip both coins and if the first coin is 0, output m as the 
outcomes of both coin tosses, but if the first coin is 1, output 
m = (1,0). Therefore there may be statistics that an IQP 
circuit would find it hard to produce that a classical computer 
could produce easily. However, it may still be possible to 
produce the above distribution using ancilla qubits (and then 
tracing them out). 



